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py, . Abstract 

Mh . In this paper, we review the theory of time space-harmonic polynomi- 

r^ ' als developed by using a symbolic device known in the literature as the 

classical umbral calculus. The advantage of this symbolic tool is twofold. 
First a moment representation is allowed for a wide class of polynomial 
stochastic involving the Levy processes in respect to which they are mar- 
tingales. This representation includes some well-known examples such as 
Hermite polynomials in connection with Brownian motion. As a conse- 
^ I quence, characterizations of many other families of polynomials having the 

^|h , time space-harmonic property can be recovered via the symbolic moment 

^p^ ' representation. New relations with Kailath-Segall polynomials are stated. 

^Nj , Secondly the generalization to the multivariable framework is straightfor- 

^Z^ • ward. Connections with cumulants and Bell polynomials are highlighted 

—^ ' both in the univariate case and in the multivariate one. Open problems 

f^ , are addressed at the end of the paper. 
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iJ '. 1 Introduction 

In mathematical finance, a Levy process |22] is usually employed to model option 
pricing. 

Definition 1.1. A Levy process X — {Xt}t>a 'is a stochastic process satisfying 
the following properties: 

a) X has independent and .stationary increments; 

b) P[X(0) = 0] = 1 on the probability space (51,7^, P); 

c) X is stochastically continuous, i.e. for alia > Q and for all s > 0, linif-s.^ P{\X{t)- 

X{s)\ >a)^Q. 



The employment of Levy processes in mathematical finance is essentially 
due to the property of manage continuous processes interspersed with jump 
discontinuities of random size and at random times, well fitting the main dy- 
namics of a market. In order to include the risk neutrality, a martingale pricing 
could be applied to options. But Levy processes do not necessarily share the 
martingale property unless they are centred. Instead of focusing the atten- 
tion on the expectation, a different approach consists in resorting a family of 
stochastic processes, called polynomial processes and introduced very recently 
in [4] . These processes are built by considering a suitable family of polynomi- 
als {P{x,t)}t>o and by replacing the indeterminate x with a stochastic process 
Xt- Introduced in [25] and called time-space harmonic polynomials (TSH), the 
polynomials {P{x,t)}t>o are such that 

(1.1) E[P{Xt, t) I 5s] - P{Xs, s), for s < t 

where S^^ = cr (^t : t < s) is the natural filtration associated with {Xt\t>Q- 

As done in [TO] for the discretized version of a Levy process, that is a ran- 
dom walk, TSH polynomials can be characterized as coefficients of the Taylor 
expansion 

E[eyi^{zXt}\ ^ fc! 

in some neighborhood of the origin. The left-hand side of ()1.2p is the so-called 
Wald's exponential martingale [1^. Wald's exponential martingale is well de- 
fined only when the process admits moment generating function £'[exp{zXt}] in 
a suitable neighborhood of the origin. Different authors have tried to overcome 
this gap by using other tools. Sengupta 25 uses a discretization procedure to 
extend the results proved by Goswami and Sengupta in [15] . Sole and Utzet [27] 
use Ito's formula showing that TSH polynomials with respect to Levy processes 
are linked to the exponential complete Bell polynomials ^3^. Wald's exponential 
martingale (|1.2p has been recently reconsidered also in [26j, but without this 
giving rise to a closed expression for these polynomials. 

The employment of the classical umbral calculus turns out to be crucial in 
dealing with (|1.2p . Indeed, the expectation of the polynomial processes Rk{Xt, t) 
can be considered without taking into account any question involving the conver- 
gence of the right hand side of (11.21) . Indeed the family {Rk{x,ty\t>o is linked 
to the Bell polynomials which arc one of the building blocks of the symbolic 
method. The main point here is that any TSH polynomial could be expressed 
as a linear combination of the family {Rk{x^t)} and the symbolic representa- 
tion of these coefficients is particularly suited to be implemented in any symbolic 
software. The symbolic approach highlights the role played by Levy processes 
with regard to which the property (jl.ip holds and makes clear the dependence 
of this representation on their cumulants. 

The paper is organized as follows. Section 2 is provided for readers unaware 
of the classical umbral calculus. We have chosen to recall terminology, notation 



and the basic definitions strictly necessary to deal with the object of this paper. 
We skip any proof. The reader interested in is referred to [lOj [11] . The theory 
of TSH polynomials is resumed in Section 3 together with the symbolic repre- 
sentation of Levy processes closely related to their infinite divisible property. 
Umbral expressions of many classical families of polynomials as TSH polynomi- 
als with respect to suitable Levy processes are outlined. The generalization to 
the multivariable framework is given in Section 4. This setting allows us to deal 
with multivariate Hermite, Euler and Bernoulli polynomials as well as with the 
class of multivariate Levy-Sheffer systems introduced in [9] . Open problems are 
addressed at the end of the paper. 



2 The classical umbral calculus. 

Let M.[x] be the ring of polynomials with real coefficientio in the indeterminate 
X. The classical umbral calculus is a syntax consisting in an alphabet A = 
{a, f3^j, . . .} of elements, called umbrae^ and a linear functional E : R[a;] [A] — > 
M.[x], called evaluation, such that E[l] — 1 and 

£;[a;" a' (3^ ■■■ 7'=] = x" E[a'] E[l3i] ■ ■ ■ £^[7'^] (uncorrelation property) 

where a, /3, . . . , 7 are distinct umbrae and n, i,j,...,k are nonnegative integers. 
A sequence {ai}i>o G M.[x], with oq = 1, is umhrally represented^ by an 
umbra a if i?[a'] = a^, for all nonnegative integers i. Then a^ is called the 
i-th moment of a. An umbra is scalar if its moments are elements of R while 
it is polynomial if its moments are polynomials of M\x\. Special scalar umbrae 
are given in Table 1. The core of this moment symbolic calculus consists in 



Umbrae 


Moments 


Augmentation e 


-B[e*] = (5o,i, with 5i,j = 1 if i = j, otherwise Si^j = 0. 


Unity u 


E[u'] = 1 


Boolean unity u 


E[u'] = i\ 


Singleton x 


E\x\ = 1 and E[x'] = 0, for alH > 1 


Ben;3 


E[/3''] = Bi, with Bi the i-th Bell number 


Bernoulli l 


E[i'] = «8i, with «8, the i-th BernouUi number 


Euler e 


E[e^] — (Bi, with €i the i-th Euler number 



Table 1: Special scalar umbrae. The equalities on the right column refer to all 
nonnegative integer i, unless otherwise specified. 

defining the dot-product of two umbrae, whose construction is shortly recalled 
in the following. 



^The ring K[x] may be replaced by any ring in whatever number of indeterminates, as for 
example M.[x,y, . . .]. 

^Whcn no misunderstanding occurs, we use the notation {fli} instead of {ai}i>o 



First let us underline that in the alphabet A two (or more) distinct umbrae 
may represent the same sequence of moments. More formally, two umbrae a 
and 7 are said to be similar when ii^[a"] = ii^[7"'] for all nonnegative integers n, 
in symbols a = 7. Therefore, given a sequence {a„}, there are infinitely many 
distinct, and thus similar umbrae, representing the sequence. 

Denote a' + a" + ■ ■ ■ + a"' by the symbol n.a, where {a', a", . . . ,a"'} is a set 
of n uncorrelated umbrae similar to a. The symbol n.a is an example of auxiliary 
umbra. In a saturated umbral calculus, the auxiliary umbrae are managed as 
they were elements of A [21]. The umbra n.a is called the dot-product of the 
integer n and the umbra a with moments jll) : 

i 

(2.1) Qiin) = E[{n.ay] = y^(n)fc_B,,fc(ai, 02, . . . , a,-fc+i), 

fe=i 

where (n)k is the lower factorial and Bik are the exponential partial Bell poly- 
nomials [3]. 

In (|2.1|) . the polynomial qi{n) is of degree i in n. If the integer n is replaced 
by i e R, in ((2?T|) we have qi{t) = J2k=iWkBi,k{ai,a2, ■ . ■ ,ai_k+i)- Denote by 
t.a the auxiliary umbra such that E[(t.ay] = qi{t), for all nonnegative integers 
i. The umbra t.a is the dot-product of i and a. A kind of distributive property 
holds: 

(2.2) {t + s).a = t.a + s.a', s,teR 

where a' = a. In particular if in ()2.ip the integer n is replaced by — i, the 
auxiliary umbra —t.a is such that 

(2.3) -t.a + t.a' = e, 

where a' = a. Due to equivalence ()2.3p . the umbra —t.a is the invers^ umbra 
of i.a. 

Let us consider again the polynomial qi{t) and suppose to replace t by an 
umbra 7. The polynomial 9^(7) G R[a;][^] is an umbral polynomial with support 
Ifl supp (9^(7)) — {7}. The dot-product of ^ and a is the auxiliary umbra 7.0 such 
that £^[(7. a)'] — E[qi{'^)\ for all nonnegative integers i. Two umbral polynomials 
p and q are said to be umbrally equivalent if E[p] — E[q], in symbols p c^ q. 
Therefore equations (|2.1|) . with n replaced by an umbra 7, can be written as 
the equivalences 

i 

(2.4) qi{^) ~ (7. a)' ~ ^(7)fc-Bvfe(ai, ^2, ■ • ■ , ai-fe+i)- 

fe=i 



•^Since — t.o and t.a are two distinct symbols, they are considered uncorrelated, therefore 
—t.a + t.a' = —t.a + t.a = e. When no confusion occurs, we will use this last similarity instead 
of (IZ3l l. 

*Tho support supp (p) of an umbral polynomial p G R[a;][y4] is the set of all umbrae occur- 
ring in it. 



Special dot-product umbrae are the a-cumulant umbra x.a and a-partition 
umbra (3. a, that we will use later on. In particular any umbra is a partition 
umbra ^11^. This property means that if {ai} is a sequence umbrally represented 
by an umbra a, then there exists a sequence {hi} umbrally represented by an 
umbra «;„, such that a = (3.k^. The umbra «„ is similar to the a-cumulant 
umbra, that is k„ = x-c^i a-nd its moments share the well-known properties of 
cumulants. |j 

Dot-products can be nested. For example, moments of (a.<r).7 can be re- 
cursively computed by applying two times formula (12.41) . Parenthesis can be 
avoided since (a.<r).7 = «.(<;. 7). In particular a./3.j, with /3 the Bell umbra, is 
the so-called composition umbra, with moments 

i 

(2.5) E[{a.p.-fy] ^^akBi^k{gi,g2,---,9i-k+i), 

fc=i 

where {ai} are moments of a and {gi} are moments of 7. When the umbra a 
is replaced by t € R, then equation (|2.5p gives the i-th moment of a compound 
Poisson random variable (r.v.) of parameter t : 

i 

Eiit.p.^Y] = J2 t''Br,k{9i,g2, ■ ■ . ,5.-fc+i)- 
fe=l 

There are more auxiliary umbrae that will employed in the following. For ex- 
ample, if E[a\ ^ 0, the compositional inverse a^^^^ of an umbra a is such that 
a.l3.a'^^^^ = a"^^^^ .f3.a = x- The derivative of an umbra a is the umbra ao 
whose moments are i?[a^] = ioi-i for all nonnegative integers i > 1. The dis- 
joint sum a+"f of a and 7 represents the sequence {ai + gi}- Its main property 
involves the Bell umbra: 

(2.6) (3.{a+j)=/3.a + p.-f. 

2.1 Symbolic Levy processes. 

The family of auxiliary umbrae {i.ajtg/, with / C M+, is the umbral counterpart 
of a stochastic process {Xt}t£i having all moments and such that -E[-'i^t] — 
E[{t.aY] for all nonnegative integers i. This symbolic representation parallels 
the well-known infinite divisible property of a Levy process, summarized by the 
following equality in distribution 

(2.7) Xt = AX,/„ + • • . + AX,/„ 



^For cumulants {Ci{Y)} of a random variable Y, the following properties hold for all 
nonnegative integers* : (Homogeneity) Ci{aY) = a^Ci{Y) for a G R, (Semi-invariance) Ci(Y+ 
a) = a + CiiY), Ci{Y + a) = C^{Y) for i > 2, (Additivity) Ci{Yi + Y2) = Ci(Fi) + Ci{Y2), if 
Yi and Y2 are independent random variables. 



with AXj/„ a r.v. corresponding to the increment of the process over an inter- 
val of amphtude t/n. The n-fold convolution p.7p is usually expressed by the 
product of n times a characteristic function E[e^^^*\ = £'[e'^'^'''-*/i]" with i the 
imaginary unit. More generally one has 



(2.8) 



E[e 



izXt] _ 



E\e 



\zXi]t 



Equation (j2.8p allows us to show that the auxiliary umbra t.a is the symbolic 
version of Xt . To this aim we recall that the formal power series 



(2.9) 



2'' 



is the generating function of an umbra a, umbrally representing the sequence 
{flj}. Table 2 shows generating functions for some special auxiliary umbrae 
introduced in the previous section. 



Umbrae 


Generating functions 


Augmentation e 


f{e,z) = l 


Unity u 


f{u,z) = e^ 


Boolean unity u 


/(S,^) = T3l 


Singleton x 


f{x.z) = l + z 


BeWp 


/(/3,z)-exp[e--l] 


Bernoulli l 


/(.,z) = z/(e--l) 


Euler rj 


/(r7,^) = 2eV[e^ + l] 


dot-product n.a 


/(n.a,z) = /(a,z)" 


dot-product t.a 


/(i.a,z) = /(a,z)* 


dot-product 7. a 


fh-a,z) = /(7,log/(a,z)) 


a-cumulant x-a 


/(x.a,^) = l + log[/(a,z)] 


a-partition p.a 


f{(3.a,z) = exp[f{a,z)-l] 


composition a./3."/ 


/(a./3.7,z) = /[a,/(7,z)- 


1] 


a-partition t./3.7 


/(t./?.7,^) = exp[t(/(7,z)- 


1)] 


derivative ao 


fiao,z) = l + zf{a,z) 



Table 2: Generating functions for some special auxiliary umbrae. 



As for infinitely divisible stochastic processes (j2.7p . one has f(t.a,z) — 
f{a,zY. It is well-known that the class of infinitely divisible distributions co- 
incides with the class of limit distributions of compound Poisson distributions 
[13] . By the symbolic method, any Levy process is of compound Poisson type [5] . 
This result is a direct consequence of the Levy-Khintchine formula |22j involving 
the moment generating function of a Levy process. Indeed, if (t>{z,t) denotes 
the moment generating function of Xt and 4>{z) denotes the moment generat- 



ing function of Xi then 4>{z^t) 



from 



From the Levy-Khintchine 



formula (f>{z) = exp[g{z)], with 

(2.10) giz) = zm+l-s^z^+ / (e^^ - 1 - zx l{|,|<i}) dKa:)). 

The term (m, s^, z^) is the Levy triplet and v is the Levy measure. The function 
4>{z,t) shares the same exponential form of the moment generating function 
/(i./3.7, z) of a compound Poisson process, see Table 2. If v admits all moments 
and if Co = m + Jr, |^. , x d(i'{x)), then the function g{z) given in ()2.10p has the 
form 

(2.11) g(^) = coz + 1 ,2^2 + /(e^--l-zx)d(Kx)). 



Thanks to (|2.1ip , the symbolic representation t.f3.j oia Levy process is such that 
the umbra 7 can be further decomposed in a suitable disjoint sum of umbrae. 
Indeed, assume 

i) (; an umbra with generating function /(<;, z) ~ 1 + z^/2, 

ii) 77 an umbra with generating function /(t;, z) = L (e^^ — 1 — z x) d(z^(a;)). 

Then a Levy process is umbrally represented by the family 
(2.12) {t./3.(coX+s?-i-77)}t>o or {i./3.(cox+sO + i-^-'7}t>o, 

due to (12.61) . Symbolic representation (|2.12p is in agreement with Ito representa- 
tion Xt — Wt+Mt+cot of a Levy process with Wt+cot a Wiener process and Mt 
a compensated sum of jumps of a Poisson process involving the Levy measure. 
Indeed the Gaussian component is represented by the symbol i./3.(cox+s<?) as 
stated in [7] , with cq corresponding to the mean and s^ corresponding to the vari- 
ance. The Poisson component is represented by i./3.r/, that is t./3.ri is the umbral 
counterpart of a random sum Sn = Yi + ■ ■ ■ + Yn, with {Yi} independent and 
identically distributed r.v.'s corresponding to 77, associated to the Levy measure, 
and N a Poisson r.v. of parameter t. The representation {t./3 .{cox+s<i+'fj)}t>o 
shows that the Levy process is itself a compound Poisson process with {Yi} 
corresponding to the disjoint sum (cox+si^-i-77). 

More insights may be added on the role played by the umbra cox+st+t?- 
Indeed the moment generating function of a Levy process can be written as 
(f>{z,t) = cxp[tg{z)] with g{z) — log (l){z). So the function g(z) in (I2.11|) is 
the cumulant generating function of Xi and 7 = cox+s'^+'H is the symbolic 
representation of a r.v. whose moments are cumulants of Xi. Therefore, in 
the symbolic representation t.a of a Levy process, introduced at the beginning 
of this section, the umbra a is the partition umbra of the cumulant umbra 
7 = coX+S'^+'7 that is a = /S.j. 

This remark suggests the way to construct the boolean and the free version 
of a Levy process by using the boolean and the free cumulant umbra |12j . 

Boolean Levy process. Let M{z) be the ordinary generating 
fmiction of a r.v. X, that is M(z) = 1+^^^-^ a^z*, where a^ = i?[X*]. 



The boolean cumulants of X are the coefficients bi of the power 
series B{z) — J2i>i ^i^^ such that M{z) = 1/[1 — B{z)]. Denote by 
a the umbra such that E[a^] = ilui for all nonnegative integers i. 
Then the umbra (/?„ such that a = u.f3.ip^ represents the sequence 
{bi} and is the a-boolean cumulant umbra. Therefore the symbolic 
representation of a boolean Levy process is t.u.j3.ip^. 

Free Levy process. The noncrossing (or free) cumulants of X are 
the coefficients r^ of the ordinary power series R{z) = 1 + X]i>i ^i^^ 
such that M(z) = R[zM{z)\. If a is the umbra with generating 
function M{z), then the a -free cumulant umbra ^s. represents the 
sequence {i\ri\. Assuming a the umbral counterpart of the incre- 
ment of a Levy process over the interval [0, 1], then the symbolic 
representation of a free Levy process is i..^„./3.(— 1.^„)^"^^. 

Some more remarks on the parameters cq and s may be added. The Levy 
process in (|2.12p is a martingale if and only if cq = 0, see Theorem 5.2.1 in [T]. 
When this happens, S[^t] = for all t > and the Levy process is said to be 
centered. Since the parameter cq allows the contribution of the singleton umbra 
X in (|2.12p , such an umbra plays a central role in the martingale property of a 
Levy process. If cq — 0, no contribution is given by x which indeed does not 
admit a probabilistic counterpart. 

If s = 0, the corresponding Levy process is a subordinator, with almost 
sure non-decreasing paths. The subordinator processes are usually employed to 
scale the time of a Levy process. This device is useful to widen or to close the 
jumps of the paths in market dynamics. Denote by Tt the subordinator process 
of Xt chosen independent of Xt. The process XT-t is of Levy type too. The 
symbolic representation of Tf is t./3.(coX+'7') so that t.l3.{cox+'r]')-l3-icox+si;+v) 
represents Xxt with r/ and r/ similar and uncorrelated umbrae. Despite its nested 
representation, the following result is immediately recovered: the process Xt^ is 
a compound Poisson process Sn with Yi a randomized compound Poisson r.v. 
of random parameter 77', shifted of cq in its mean. 



3 Time-space harmonic polynomials. 

Set X — {a}. The conditional evaluation E{- \ a) with respect to a handles the 
umbra a as it was an indeterminate ,8 . In particular, E{- \ a) : M.[x][A] — > M.[X] 
is such that £^(1 | a) = 1 and 

S(a;™a"7*^^' • • • | a) = x"'a''E[-/']E[e] ■ ■ ■ 

for uncorrelated umbrae a, 7, ^, . . . and for nonnegative integers to, n, i^ j, . . . . 
As it happens in probability theory, the conditional evaluation is an element of 
K.[a;][yl] and, if we take the overall evaluation of E{p \ a), this gives E[p], with 
p G M.[x][A], that is E[E{p \ a)] = E[p]. Umbral polynomials p, not having a in 
its support, are such that E{p \ a) ~ E[p]. 



Conditional evaluations with respect to auxiliary umbrae need to be handled 
carefully. For example, since (n + l).a = n.a + a' , the conditional evaluation 
with respect to the dot product n.a is defined as 

E[{n + l).a I n.a] — n.a + E[a'], 

and more general, from (j2.2p with t and s replaced by n and m, 
(3.1) 

E{[{n + m).af \ n.a) = E{[n.a + m.a'f \ n.a) = ^ H ) {n.a)^ E[{m.a')''-^], 

for all nonnegative integers n and m. Therefore, for i > the conditional eval- 
uation of t.a with respect to the auxiliary umbra s.a, with < s < i, is defined 
according to (I3.ip such as 

Elit.a)'' I s.a] = J2 (^) i^-ayEiKt - s).af-^). 
j=o ^-^^ 

Equation (jl.ip traces the way to extend the definition of polynomial processes 
to umbral polynomials. 

Definition 3.1. Let {P{x,t)} G M[x][,4] be a family of polynomials indexed 
by t > 0. P{x,t) is said to be a TSH polynomial with respect to the family of 
umbral polynomials {q(t)}t>o if and only if E [P{q{t), t) \ q{s)] = P{q{s), s) for 
alio < s <t. 

The main result of this section is the following theorem [8] . 
Theorem 3.2. For all nonnegative integers k, the family of polynomial]^ 

Qkix,t)^E[{x-t.a)''] €R[x] 
is TSH with respect to {t.a}t>o. 

By expanding Qk{x,t) via the binomial theorem, one has 

j=o ^-^^ 
so that 

Qk{t.a,t) = J2[^^it-^yE[{-t.a)''-^]. 
Since t.a is the symbolic version of a Levy process, the property 

E[Qk{t.a,t) I s.a] = ^ r/\E[{t.ay\s.a]E[{~t.a)''-^] ^Qk{s.a,s) 



^When no confusion occurs, we will use the notation x — t.a to denote the polynomial 
umbra —t.a + x = x + {—t).a. 



10 



parallels equation (|l.ll) . In particular {Qk{x,t)} is a polynomial sequence um- 
brally represented by the polynomial umbra x — t.a, which is indeed the TSH 
polynomial umbra with respect to t.a. Polynomial umbrae of type x + a are 
Appcll umbrae [6]. Then {Qk(x,t)} is an Appell sequence and 

—— Qk{x,t) — kQk-i{x,t), for all integers fc > 1. 

da; 

The generating function of the TSH polynomial umbra x — t.a is 

•> ^ ' ' k>0 

By replacing x with t.a in (|3.2I) . Wald's exponential martingale (|1.2p is recov- 
ered. Equality of two formal power series is given in terms of equality of their 
corresponding coefficients, so that E[Rk{Xt,t)] — E[Qk{t.a,t)] by comparing 
dXH) with ([12]). Wald's identity J2k>o ^i^k{Xt,t)]z'' /k\ = 1 is encoded by the 
equivalence t.a — t.a = e obtained from x ~ t.a when x is replaced by t.a. 

The next proposition gives the way to compute the coefficients of Qk{x,t) 
in any symbolic software. 

Proposition 3.3. // {a„} is the sequence umbrally represented by the umbra 
a and {Qkix,t)} is the sequence of TSH polynomials with respect to {t.a}t>o, 
then 

k 



„w,.^. 



with 






J2 d,(-l)2'W+^s[^(A),^]a^^a^^ 

Ahfc-j 



where the sum is over all partitionJij A = (f^ , 2^^ , . . .) h fc — j, s[?(A), i] denotes 
a Stirling number of first kind and d\ — i!/(ri!r2! ■ • • (l!)''^ (2!)''^ ' ' ' )• 

More properties on the coefficients of Qk{x,t) are given in [8]. 

Any TSH polynomial is a linear combination of {Qkix, t)}, which indeed are 
a bases of the space of TSH polynomials. The following theorem characterizes 
the coefficients of any TSH polynomial P{x,t) in terms of the coefficients of 
{Qk{x.t)}. 



^Recall that a partition of an integer j is a sequence A = (Ai, A2, . . . , Am), where Aj are 
weakly decreasing positive integers such that JZTLi -^j = *■ The integers \j are named parts of 
A. The length of A is the number of its parts and will be indicated by /(A). A different notation 
is A = (fi , 2'"2 ,...), where Tj is the number of parts of A equal to j and ri + r2 + ■ ■ ■ = '(A). 

o be tl 
"A is a partition of i" 



11 



Theorem 3.4. A polynomial P{x,t) — 'Ylj^oPji^)^'' "/ degree k for all t > 
is a TSH polynomial with respect to {t.a}t>Q if and only if 

pjit) ^ E (') Prwm-t-c^y-'], for J = 0, . . . ,fc. 

3.1 Cumulants. 

A different symbolic representation of TSH polynomials {Qk{x,t)} is 

Qk{x,t)=E[{x-t.l3.Kj''], 

with Kc the a-cumulant umbra. The umbra —t.p.n^ = i./3.(— I.Kq) = i./3.K(„i „) 
is the symbolic version of a Levy process with sequence of cumulants of Xi 
umbrally represented by K(_i^). Therefore, also the polynomials Qk{x,t) — 
E[{x + t.p.K^)^] are TSH with respect to Levy processes umbrally represented 
by {i./3.«;(_i.„)}t>o = {t.{-l.a)}t>o. 

Moments of polynomial umbrae i./3.7 involve the exponential Bell polyno- 
mials. When t is set equal to 1, then complete exponential Bell polynomials are 
recovered. More generally, the moments of a: — i./J.^^ = x — t.a can be expressed 
by using exponential complete Bell polynomials too since 
(3.3) 

X-t.a = [i\X-{x - t.a)] = /?.K,_t.„ = /3.(K(a;.«)+K(-t.a)) = P-K{x.u) + P-K.(-t.a) 

where K^^t.a is the cumulant umbra of x — t.a, that could be replaced by 
K{x.u)+K(-t.a) due to the additivity property of cumulants. The last equiva- 
lence in p.3p follows from equivalence (|2.6p . From equivalences p.3p . we have 

(3.4) Qk{x, t) = Yk{x + hi,h2, ..., hk), 

with Yk exponential complete Bell polynomials and {hi} cumulants of —t.a. 
Equation (|3.4p has been proved in [27] by using Teugel martingales. 
For Qk{x,t), the ShefFer identity with respect to t holds: 

Qk{x,t + s) = "^( .jPj{s)Qk-j{x,t), 
j=o ^^^ 

where Pj{s) — Qj{0, s) for all nonnegative integers j. 

3.2 Examples. 

The discretized version of a Levy process is a random walk Sn = Xi + X2 + 
• • • -|- Xn, with {Xi} independent and identically distributed r.v.'s. For the 
symbolic representation of a Levy process we have dealt with, the symbolic 
counterpart of a random walk is the auxiliary umbra n.a. Indeed the infinite 
divisible property p.7p is highlighted in the summation a' + a" + ■ ■ ■ + a'" , 
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encoded in the symbol n.a, with a' , a" , . . 



.,a 



uncorrelated and similar umbrae. 



Nevertheless not all r.v.'s having the symbolic representation n.a share the 
infinite divisible property. For example, the binomial r.v. has not the infinite 
divisible property [22j . nevertheless its symbolic representation is of type n.a 
where a = X-P-l^ andp G (0, 1). So the generality of the symbolic approach lies in 
the circumstance that if the parameter n is replaced by i, that is if the random 
walk is replaced by a Levy process, more general classes of polynomials can 
be recovered for which many of the properties here introduced still hold. The 
following tables resume the TSH representation for different families of classical 
polynomials, see [20]. In particular Table 3 gives the umbra corresponding to 
the r.v. Xi of 5„ in the first column, its umbral counterpart in the second 
column and the associated TSH polynomial in the third column. In Table 4, 
the TSH polynomials given in Table 3 are traced back to special families of 
polynomials. In particular, with the polynomials Vkixjt) we refer to 

k 
Vk{x,t) = y^ Qj{x, t)Bkj (toi, 7712, ■ ■ ■ , nik-j + l) 

for suitable {Qj{x,t)} and {nii]. 



X, 


Umbral counterpart 


Corresponding TSH polynomial 


Uniform [0, 1] 


-l.i 


E[{x + n.Lf] 


Bernoulh p= \ 


\{-l.e + u) 


E[{x + n.[\{-l.u + f)\f] 


Bcrnoulh p G (0, 1) 


X-P-P 


E[{x - n.x-p.fif] 


Sum of a e N 
uniform r.v.'s on [0, 1] 


a.{-l.L) 


E[{x + {an).Lf] 



Table 3: TSH polynomials associated to special random walks 



X, 


Special families 
of polynomials 


Connection with 
TSH polynomials 


Uniform [0, 1] 


Bernoulli Bk{x,n) 


Bk{x,t) = Qk{x,t) 


Bernoulli p = 1/2 


Euler £k{x^ n) 


£k{x,n) = Qk{x,t) 


Bernoulh p£ (0, 1) 


Krawtchouk ]Ck{x,p, n) 


{n)klCkix,p,n) = Pkix,t) 
m,=E[ii-l.x.p.l3)<-'>y] 


Sum of a e N 
uniform r.v.'s on [0, 1] 


pseudo-Narumi Nk{x, an) 


k\Nk{x,an) ^Vk{x,t) 
m, ^ E[{u<-^>Y] 



Table 4: Connection between special families of polynomials and TSH polyno- 
mials 



Next tables 5 and 6 give TSH polynomials for some special Levy processes. 
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Levy process 


Umbral representation 


TSH polynomial Qk{x,t) 


Brownian motion 
with variance s^ 


t.(i.{s^) 


E[ix-t.(3.{.s,m 


Poisson process 
with parameter A 


(U)./3 


E[{x {t\).m 


Gamma process 

with scale parameter 1 

and shape parameter 1 


t.u 


E[{x-t.u)''] 


Gamma process 

with scale parameter A 

and shape parameter 1 


{t\).u 


E[{x - itX).u)''] 


Pascal process 

with parameter d — p/q 

and p + q = 1 


t.u.d.fi 


E[{x-t.u.d.l3f] 



Table 5: TSH polynomials associated to special Levy processes 



3.3 Orthogonality of TSH polynomials. 

A special class of TSH polynomials is the one including the Levy-ShefFer poly- 
nomials, whose applications within orthogonal polynomials are given in ^3]. A 
sequence of polynomials {Vk{x, t)}t>Q [M] is a Levy-ShefFer system if its gener- 
ating function is such that 



(3.5) 



'^yk{x,t)— = (.g(z))*exp{xu(z)}, 



fc>0 



where g{z) and u{z) are analytic functions in a neighborhood of z = 0, m(0) = 0, 
g(0) = 1, m'(0) ^ and l/g{T{z)) is an infinitely divisible moment generat- 
ing function, with t(z) such that t{u{z)) — z. Assume a an umbra such that 
f{a,z) = g{z) and 7 an umbra such that f{'^,z) = 1-1- u{z). From p.Sp . the 
Levy-ShefFer polynomials are moments of x./3.^ + t.a : 



(3.6) 



Vkix,t)=E[{x.p.j + t.a)'']. 



Theorem 3.5. The TSH polynomials Qk{x,t) are special Levy-Sheffer polyno- 
mials. 



The proof of Theorem 13.51 is straightforward by choosing in p.6p as umbra 
a its inverse —l.a and as umbra 7 the singleton umbra x- AH the Levy-Sheffer 
polynomials possess the TSH property. Indeed the following theorem has been 
proved in [8]. 

Theorem 3.6. The Levy-Sheffer polynomials {Vk{x,t)}t>Q are TSH with re- 
spect to Levy processes umhrally represented by {— t.a./3.7'^^^''}f>o. 
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Levy process 


Special family 
of polynomials 


Connection with 
TSH polynomials 


Brownian motion 
with variance s^ 


HermiteiJ^''^(a;) 


Hf\x) = Qk{x,t) 


Poisson process 
with parameter A 


Poisson-Charlier Ck{x, Xt) 


Ckix,Xt) = 
= T!;^is{k,j)Qu{x,t) 
with s{k,j) Stirling 
numbers of first kind 


Gamma process 

with scale parameter 1 

and shape parameter 1 


Laguerre Cl^''{x) 


kl{-l)''Cl''{x)^Qk{x,t) 


Gamma process 

with scale parameter A 

and shape parameter 1 


actuarial gk{x,t) 


gkix,t)^Vkix,t) 


Pascal process 

with parameter d = p/q 

and p + q = I 


Meixner polynomials 
of first kind Mi^{x,t,p) 


{-lf{t)kMk{x,t,p)^ 
^Vk(x,t) 
mi=E[{x.{^l.X + x/pm 



Table 6: Special families of polynomials and TSH polynomials 



In particular, one has [8] 

k 

(3.7) Vu{x, i) - XI ^[("^ + ^■/3-'^(a./3.7<-i>))'i^fc,»(ffi' ■ • • ' 9k-^+l). 



i=0 



where gi = £'[7*], for all nonnegative i, and K(Q./3.'y<-i>) is the cumulant umbra 
of a./3.7^^^^, with ^<"^> the compositional inverse of the umbra 7. When the 
umbra a is replaced by its inverse and the umbra 7 by the singleton umbra, since 
X^"^> = X, the only contribution in the summation p.7p is given by « = fc. So 
again equation (13.71) reduces to Qk{x,t) = E[{x — t.a)'^] since K(a.^.7<-i>) = 
X-- 1-a- 

Within Levy-Sheffer polynomials, the Levy- Meixner polynomials are those 
orthogonal with respect to the Levy processes —t.a.fS.'^'^^'-^, due to their TSH 
property. The orthogonal property is 

According to [23] , all the polynomials in Table 6 are orthogonal. Their measure 
of orthogonality corresponds to the Levy process —t.a since x^^^^ = X ^-nd 
-t.a.l3.j<-^> = -t.a. 
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3.4 Kailath-Segall polynomials. 

Equivalence p.3|) gives the connection between TSH polynomials and Kailath- 
Segall polynomials (16j , which is a different class of polynomials strictly related 
to Levy processes. Indeed, both have a representation in terms of partition um- 
bra of a suitable polynomial umbra. Overlaps are removed by suitably choosing 
the indeterminates. 

The n-th Kailath-Segall polynomial Pn{xi, . . . , x„) is a multivariable polyno- 
mial such that when the indeterminates are replaced by the sequence X^ , . . . , 
Xf of variations of a Levy process Xt 

X[^^ = Xu Xf ) - [X,X]u Xt^ = E(^^«)" " ^ 3, 

s>t 

its iterated integrals are recovered 

pW = l, P,^'^=Xt, p/"'= Tpir'^dX,, n>2, 



that is P( — Pn [Xf. , . . . , Xf J . The following recursion formula is known 
as Kailath-Segall formula 

(3.8) P/") = - (P/""'^X,^') - p/"-'^xf ^ + • • . + (-l)"+ip/"^X,(")) . 

When Xf , . . . , X^ are replaced by the power sums Si, . . . , Sn in the indeter- 
minates xi, . . . ,Xk, according to formula (1.2) in [28] and Theorem 3.1 in [13], 
the corresponding polynomials Pn(S'i, . . . , 5'„) are such that 

(3.9) n!P„(5i, . . . , 5„) = P[(/3.[(x.x)fT])"], 

where a is the power sum umbra representing {Sj} and the x-cumulant umbra 
X-X represents the sequence {(— l)*^^(i — 1)!}. 

In order to recognize special TSH polynomials within the family {Pn}, two 
steps are necessary: 

i) Kailath-Segall polynomials need to be umbrally represented when 
the power sums {Sj} are replaced by the indeterminates {xi}; 
ii) the indeterminates {xi} need to be replaced by suitable terms 
involving x and t. 

For the first step, we will use equation p.9p . Assume p an umbra representing 
the sequence {E[{x-Xi.(3y]}. Observe that E\{x-Xi.l3y] — xi for all nonnegative 
i. Then from p.9p one has n!Pn(xi, . . . ,a;„) — i?[(/3.[(x.x)p])"] so that the 
generating function of P„ is 

/(/3-[(x-x)p],^)=exp j^i^^z"a;„J , 
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see also [30]. The strength of this symboUc representation essentially relies 
on the properties of the partition umbra [i\{x-x)p] reproducing those of Bell 
polynomials. For example, the following property of Kailath-Segall polynomials 

(3.10) F„(aa;i,a^X2,. . . ,a"a;„) = a"F„(a;i,X2, . . . ,x„), aGM 

is proved by observing that j3.[a(x-x)p] = a(/5-[(x-x)p])- For the next step, we 

need to characterize the indeterminates xi^X2, ■ ■ ■ such that 

(3.11) 

^i{x.u) + ^(-t.a) = {X-X)P ^ ^[('«(:..u))"] + ^[(K(_t.a))"] = (-1)" ^{n - 1)! X„. 

In the following we show some examples of how to perform this selection. These 
results extend the connections between TSH polynomials and Kailath-Segall 
polynomials analyzed in [27 . 

Generalized Hermite polynomials: 

Since ^[(k(^.„))*] = xJj,i and ^[(K(_t.;3. («?)))'] = s^i<5i,2, then 

k\Pk{x,s\Q,...,Q)=H'^\x) 

where Y.k>n Hi^\x)z''/kl = exp{xz - tz^/2}. 
Poisson-Charlier polynomials: 

Poisson-Charlier polynomials {Ck{x,t)}, with generating function 

Y.Ckix,t)^ = e''%l + zr, 

fc>0 

are umbrally represented by 

(3.12) Ckix,\t)^E[{x.x-t.\.u)''], 

see [S]. Nevertheless (|XT^ differs from the resuh of Theorem E21 the TSH 
property holds since {Ck{x, Xt)} are a linear combination of special Qk{x,t). 
Moreover representation (|3.12p allows us the connection with Kailath-Segall 
polynomials, when the indeterminates {xi] are chosen such that E[{Hi(^,j. ^-^y] + 
E[{n(^t.x.u)Y] = i-iy-H^ - 1)! a:,. Since 



f X — tX 

^[(«(x.x))1 + ^[('«(-t.A.n))1 = I (_l)«-l(j _ 1)1 



X — tX, * = 1, 

i > 2, 



then klPk{x — tX, x,x, . . .) = Ck{x, Xt). 
Laguerre polynomials: 

Laguerre polynomials {jO.*f7 (x)} arc TSH polynomials such that 

k\i-i)'ci-'{x)^E[ix-t.un 

They can be traced back to Kailath-Segall polynomials if the indeterminates 
{xi} are characterized by i?[(K(j, „))*] + -E[(k(_j g))*] = (— l)*^-^(i — 1)1 Xi. Since 
/('*(-*.«), ^) = 1 + i log(l - z) then £;[(K(_,.a))*] = -t{i - 1)!. So 

Pk{x-t,t,~t,t,...) = {-!)'' Cl'^ix) 
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and from ([5TII)) we have Pk{t ~x,t,t,...)^ ^V^i^)- 
Actuarial polynomials: 

The actuarial polynomials gk{x,t) are a linear combination of suitable TSH 
polynomials Qk{x,t) (see Table 6) but they are moments of the polynomial 
umbra Xt — x.(3, that is gk{x,t) = E[{Xt — x.(3)'^]. In order to characterize the 
connection with Kailath-Segall polynomials, the indeterminates {xi} need to 
be characterized by E[{K(^xt.u)T] + -^[(«^(-2;.u))'] == (^1)* H* ^ ^V-Xi. As be- 
fore E[(K(^\t,u)y] = Xtdi^i, instead E[{K(^-x.u)y] = —x{i — 1)! for all nonnega- 
tive integers i as in the previous example. Therefore one has fc!(— l)'^Pfe(a; — 
Ai,x,x, ...) = gk{x,t). 

Meixner polynomials of first kind: 

Meixner polynomials of first kind {Mfe(x, t,p)} [23] are a linear combination 
of suitable TSH polynomials Qk{x, t) (see Table 6) but they are moments of the 
following polynomial umbra 



i-l)''{t)kMkix,t,p)^E- 



-i-x 



t.x 



which allows us to find the connection with Kailath-Segall polynomials. Indeed 
for all nonnegative integers i we have 



{''-(-i-x+f)} 



{-ly-'{^-ly.x 



1 



and 



E 



{^(x-t.^)}' =(-ir'(*-i)!i. 



Then Kailath-Segall polynomials give the Meixner polynomials {—!)'' {t)kMk{x, t,p) 
by choosing 

— -1] x-t 
for i = 1,2,.... 



4 Symbolic multivariate Levy processes. 



In the multivariate case, the main device of the symbolic method here proposed 
relies on the employment of multi-indices of length d. Sequences like {gii,i2,...,id} 
are replaced with a product of powers fJ'i 1^2 ' ' ' Md*' where (/ii, /i2, . . . , /id) are 
umbral monomials and (ii, 12, . . . , id) are nonnegative integers. Since the umbral 
monomials could not have disjoint support, then the evaluation E does not 
necessarily factorizes on the product i^-l 11^2 '' ' Ki > ^^^^^ ^^ 



(4.1) 



E[^^^^ 



■Mdl 



Eifj,' 



9i 



where i — (ii, 12, . . . , id) and /x = (/zi,/^2, • ■ • ,Md)- We assume go — 1 with 
= (0,0,..., 0). Then gi is called the multivariate moment of /x. Table 7 shows 
some special d-tuples we will use later. 
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d- tuple 


Generating functions 


Multivariate Unity u 


{u,...,u') 


fiu,z) =6^1+-+^''. 


Multivariate Gaussian <; 


(C,...,0 


f{<;,z) = l + ^^zz^. 


Multivariate Bernoulli t 


(t, . . . , t) 




Multivariate Euler q 


iv.---.r]) 





Table 7: Generating functions of special d-tuples of umbral monomials 



The notions of similarity and uncorrelation are updated as follows. Two 
d-tuples /x and v of umbral monomials are said to be similar if they represent 
the same sequence of multivariate moments. They are said to be uncorrelated 

Multivariate Levy processes are represented by d-tuples of umbral monomi- 
als. 



Definition 4.1. A stochastic process {Xt}t>o on 
process if 



is a multivariate Levy 



i) Xo = a.s. 

ii) For all n > 1 and for all < ii < ^2 < 
Xti , Xt3 — Xt2 , ■ ■ ■ are independent. 



< tn < oo, the r.v.'s Xt 



iii) For all s < t, Xt-i^s ^ Xg = Xt- 

iv) For all e > 0, lim,i^o PQXt+h - Xt| > e) := 0. 

v) i I— > Xt{uj) are right- continuous with left limits, for all w G i7, with CI the 
underlying sample space. 

As in the univariate case, the moment generating function of a multivariate 



Levy process is (pxj^{z) = E 



with 2 6 R''. Paralleling the univariate 



case, the generating function of a d-tuple /x is 



/(M^^) = i+EE54 



fe>i , 



\i\=k 



Choose the d-tuple /x such that f{ij,,z) — (pxi{z), that is E[fi'^] — E[X\] for 
all i E Nq. The auxiliary umbra n.fi denotes the sum of n uncorrelated d-tuples 
of umbral monomials similar to /x. Its multivariate moment is [5] 



(4.2) 



E[in.^,y] ^ Y. 



Ahi 



a(A)A! 



{n)i^x) E[nx]: 



19 



where E[fi^] ~ g^ g^ . . . and A is a partitioro of the multi-index i of length /(A). 
By replacing the nonnegative integer n with the real parameter t in (|4.2I) the 
resulting auxiliary umbra t.fi is the symbolic representation of the multivariate 
Levy process Xt- 

As in the univariate case, since fi = l3.Kfj, with k^ the /x-cumulant umbra 
|13| . a different representation for a multivariate Levy process is i./3.K^. The 
cumulant d-tuple could be further specified by using the multivariate Levy- 
Khintchine formula [22|. 

Theorem 4.2. X ~ {Xt}t>o is a Levy process if and only if there exists 
a symmetric, positive defined dx d matrix E > and a measure v on 



mi e 



with 



such that 
(4.3) 

^xiz) =exp<i t 



i/({0}) = and / {\x\^ A l)iy{dx) < oo 



— zT,z^ + miz^ 



(e" 



1 - a;2;'^l{|^|<i}(a;)) i^(da;) 



The representation of ipx{z) in (RTM) by rrii, T, and v is unique. 
Set m2z'^ — /jgrf 2;a;'^l{|a,|>ii,(a;) ^{dx) and m = mi + m2, then 
'1 



y^xiz) =exp<^i 
that is, 
(4.4) ip^{z) =expii 



'-zJ:z 



T , T 

' mz 



—z'Ez'^ +mz'^ 



[e""^ -l-zx')v{dx) 



cxp -^ i / (e^^ -\-zx'^)i'{dx) 



Theorem 4.3. Every Levy process {Xt}t>o on M'' is umhrally represented by 
the family of auxiliary umbrae 



(4.5) 



{t./3.(x.m+^C^+77)} 



t>Q: 



where /3 is the Bell umbra, m G M"^', ^ is the multivariate umbral counterpart of 
a standard gaussian r.v., C is the square root of the covariance matrix S and r) 
is the multivariate umbra associated to the Levy measure. 

Every auxiliary umbra t.fi.n^ is the symbolic version of a multivariate com- 
pound Poisson r.v. of parameter t, that is a random sum Sn = Yi + - ■ ■ + Y" jv of 
independent and identically distributed random vectors {^i}, whose index N 



*A partition A of a multi-index i, in symbols A h i, is a matrix A = (A^j) of nonnegative 
integers and with no zero columns in lexicographic order -i such that A^ +A,.2 +■ ■ ■ + ^r,. = ir 
for r = 1,2, ... ,d. The number of columns of A is denoted by /(A). The notation A = 
(Aj^, Aj^, . . .) represents the matrix A with r\ columns equal to Ai, r2 columns equal to 
A2 and so on, where Ai ^ A2 -< ••.• We set m(A) = (ri,r2,...), m(A)! = ri!r2!-- and 
A! = Ai!A2!--- 
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is a Poisson r.v. of parameter t. Then the same holds for the Levy process. The 
d-tuple (^C-^-i-x.m-i-r/) umbrally represents any of the random vectors {^i}- 
Observe that x.m has not a probabihstic counterpart. If m is not equal to 
the zero vector, this parallels the well-known difficulty to interpret the Levy 
measure as a probability measure. 

4.1 Multivariate TSH polynomials. 

The conditional evaluation with respect to an umbral d-tuple /x has been in- 
troduced in [9;. Assume X — {/ii, fi2, ■ ■ ■ , fJ-d}- The conditional evaluation with 
respect to the umbral d-tuple /x is the linear operator 

E{-\fi): R[xi,...,Xd][A] -^ R[X] 

such that E{1 \ fj.) = 1 and 

(4.6) E{x[' x^i ■■■ x^l /x^ v^^^■■■\^l)^ x\ x'i ■■■ x^^ fi' Ey] E['n''] ■ ■ ■ 

for uncorrelated d-tuples fi,v>,r] . . . , multi-indices i,j,k... G Ng and {h}f^i 
nonnegative integers. Since f[{n+m).fi, z] = f{fi,z)^^™ = f{n.fi,z)f{m.fi,z), 
then {n -\- m).^ = n.fj, + m.fi' , with /i and /x' uncorrelated d-tuples of um- 
bral monomials. Then, for E{ ■ \ n.fi) we assume E[{{n + m).fiY I n.fi] — 
E[{n.fi + m.fi'Y I n.fi] for all nonnegative integers n, m and for all i G Ng. If 
n ^ m, then 

(4.7) E[{{n + Tn).fiy \ n.fi] = E[{n.fi + m.fiY \ n.fj.], 

since n.fj, and m.fj, are uncorrelated auxiliary umbrae. We will use the same 
d-tuple /x as in (j4.7p when no misunderstanding occurs. Thanks to equations 
(|46| and (|4J| . we have 

[ 

k<i 



(4.8) E [{{n + m).f,r \ n./x] = ^ i\{n.f,fE[{m.^y-% 



where k < i ^ kj < ij ioi all j ^ 1, . . . ,d and (^) = (i') ' ' ' (i'') ■ ^y analogy 
with (|4.7p and (|4.8|) . we have t.fi = s./x + {t — s)./i and for t > and s <t 



fc<i ^ ^ 

Theorem 4.4. For all i G Nq, t/ie family of polynomials 
(4.9) Qi(a;, ^ = E[{x - t.fiY] G R[xi, ...,Xd] 

is TSH with respect to {t.^}t>o. 



1 i—k] 
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The auxiliary umbra — i./x denotes the inverse of t.fj, that is —t.fi + t./j, = 
e, where e is the d-tuple such that e — (ei,e2, . . . , e^) with {ei} uncorrelated 
augmentation umbrae. Coefficients of Qi{x,t) in (|4.9p are such that 



k<i 



Q,(a.,i)=^(Ma.^-'=i?[M./x)1 



so when x is replaced by t.fi their overall evaluation is zero. Properties on the 
coefficients of Qi{x, t) can be found in [9 . Here we just recall a characterization 
of the coefficients of any multivariate TSH polynomial in terms of those of 
Qi{x,t). 

Theorem 4.5. A polynomial 

(4.10) F(a;,t) = ^P.Wa='" 

k<v 

is a TSH polynomial with respect to {t.fi}t>o ^/ o.nd only if 

(4.11) p,(i)= Y. (l) p.m m-t.(j^r-% fork<v. 

Table 8 and 9 give some examples of multivariate TSH polynomials and their 
connection with multivariate Levy processes. The corresponding d-tuples are 
given in Table 7. 



multivariate Levy process 


Umbral representation 


TSH polynomial Qi{x,t) 


Brownian motion 

with covariance S = CC'^ 


t.p.{<;C^) 


E[{x~t.p.{c;CT)f] 


Xt with Xi = ([/,..., t/) 
and U uniform r.v. on [0, 1] 


-t.L 


E[{x + t.LY] 


Xt withXi = (r,...,y) 

and Y BcrnouUi r.v. 
of parameter 1/2 


i[i.(«-l.r,)] 


E[{x + \Urj-u)]y] 



Table 8: TSH polynomials associated to special multivariate Levy processes 

Let us remark that Hermite polynomials Hi '{x,Yi) in Table 9 are a gen- 
eralization of the polynomials Hi{x) in [53] whose moment representation is 
Hi{x) = E[{xY.-^ +\YY] with E the expectation symbol, Y ~ iV(0, E"^) and 
S a covariance matrix of full rank d. 

A generalization of Levy-Sheffer system to the multivariate case has been 
introduced in [S|. A sequence of multivariate polynomials {Vfc(a;, t)}f>o is a 
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multivariate Levy process 


Special family 
of polynomials 


Connection with 
TSH polynomials 


Brownian motion 

with covariance E = CC'^ 


YLciuAic Hf\x,T.) 


Hf\x,^) = Qi{x,t) 


Xt with Xi^{U,...,U) 
and U uniform r.v. on [0, 1] 


BcrnoulhBf'(a;) 


Bf\x) = Qi{x,t) 


Xt with Xi = {¥,...,¥) 
and Y Bcrnouhi r.v. 
of parameter 1/2 


¥.n\ci £f\x) 


£f{x) = Qi{x,t) 



Table 9: Special families of polynomials and TSH polynomials 
multivariate Levy-ShefFer system if 

k 

1 + ^ ^ Vk{x,t)^ = [g{z)fcx:p{{xi + --- + Xd)[h(z)-l]}, 

\v\=k 

where g{z) and h{z) are analytic in a neighborhood of z = and 
d 



d: 



-h{z] 



^0 fori = 1,2, 



If fj, and u arc d-tuples of umbral monomials such that f{p,,z) = g{z) and 
f{v,z) = h{z) respectively, then 



(4.12) 



Vk{x, t) = E[{t.^l +{xi + --- + Xd).li.vf] 



The multivariate Levy-Sheffer polynomials for the pair /x and v are TSH polyno- 
mials with respect to a special symbolic multivariate Levy process involving the 
multivariate compositional inverse of a d-tuple u. Assume Xr^\ the rf-tuple with 
all components equal to the augmentation umbra and only the i-th one equal 
to the singleton umbra, that is Xu) = {^t ■ ■ tXt ■ ■ t^)- The multivariate com- 
positional inverse of u is the umbral rf-tuple u^'^^ = ((i/<"^>)j^, . . . , {v^^'^^)d) 
such that {u'^^^^)i.j3.u = X{i) for i ~ 1, . . . ,(i. 

Theorem 4.6. The multivariate Levy-Sheffer polynomials for the pair fj. and 
u are TSH polynomials with respect to the symbolic multivariate Levy process 

{t.{ii^.f3.v<-'> + • • • + /irf./?.i.<-^>)}t>o. 



5 Conclusions and open problems. 

In this paper, the review of a symbolic treatment of TSH polynomials, relied 
on the classical umbral calculus, is proposed. The main advantage of this sym- 
bolic presentation is the plainness of the overall setting which reduces to few 
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fundamental statements, but also the availability of efficient routines [7] for the 
implementation of formulae as (J4.2I) , which is the key to manage the polynomials 
Qk{x,t). 

The main result of this presentation is that any univariate (respectively mul- 
tivariate) TSH polynomial has the form Qk{x,t) (respectively Qk{x,t)) or can 
be expressed as a linear combination of the polynomials Qk{x,t) with coeffi- 
cients given by (j4.1ip . Thanks to the umbral representation of multivariate 
Levy-Sheffer systems, more families of umbral polynomials could be character- 
ized, together with their orthogonality properties. This will be the object of 
future research and investigation. 

In [2], Barrieu and Shoutens have related the infinitesimal generator of a 
Markov process to a more general class of linear operators possessing the TSH 
property, both ascribable to special families of martingales. A stochastic Taylor 
formula is produced which results to be a generalization of a TSH polynomial 
due to the presence of a remainder term series. A symbolic representation of this 
new TSH function could open the way to a new classification of the correspond- 
ing operators by which to recover the martingale property on Levy processes. 
Similarly, the extension to the more general class of Markov processes (a first 
attempt is given in [2] ) would move the employment of TSH functions beyond 
the field of applications strictly connected to the market portfolio. One step 
more consists in dealing with matrix-valued stochastic processes by replacing 
formal power series (|2.9|) with hypergeometric functions, as done in |18j . This 
would allow us a symbolic representation also for zonal polynomials whose com- 
putational handling is still an open problem. 
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